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1.  INTRODUCTION 


The  purpose  of  this  research  Is  to  determine  the  radar  scattering 
characteristics  of  an  airborne  target  with  complex  structure  such  as  an 

airplane.  Even  though  this  subject  matter  has  been  studied  extensively 

•/ 

in  the  past  twenty  years  [1,  2] (there  is  no  one  method  that  is  capable 
of  producing  efficient  and  accurate  results  under  all  conditions.  One 
has  to  devise  different  computational  schemes  and  employ  different 
mathematical  techniques  which  depend  on  factors  such  as  the  particular 
geometry  of  the  target,  the  frequency  range  of  interest,  the  acceptable 
error,  and  computation  labor. 

In  the  present  study,  we  intend  to  develop  a  relatively  simple 

analytical  model  of  scattering  from  modern  airplanes  for  the  frequency 

range  between  1  and  10  GH  .  Great  precision  of  the  theoretical  model  is 

z 

not  required  since  only  the  smoothed  curves  of  the  calculated  data 
(e.g. ,  medians  over  small  angular  region)  will  be  of  interest.  Pre¬ 
liminary  calculation  is  to  be  kept  to  a  minimum  because  repeated 
calculations  for  a  wide  range  of  parameters  will  be  performed  for 
later  statistical  studies  in  the  simulation  of  a  particular  radar 


system  [2] 
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With  the  above  objective  in  mind,  ve  have  searched  the 'literature 


and  evaluated  the  different  techniques.  It  has  been  concluded  that  the 
ray-optical  method  [3-19]offers  the  best  approach  for  the  following 

reasons: 

(i)  This  method  is  capable  of  handling  target  sizes  that  are  large 
compared  to  the  wavelength. 

v’ii)  The  ray-optical  method  yields  a  final  result  for  the  scattered 
field  in  terms  of  simple  functions,  thus  allowing  efficient 
numerical  computation. 

(iii)  The  ray-optics  approach  offers  a  convenient  means  for  taking 
the  shadowing  effect  into  account.  This  is  accomplished  by 
first  locating  the  scattering  centers,  or  the  bright  points, 
and  adding  a  simple  logic  in  the  program  to  indicate  which  of 
the  bright  points  are  in  the  shadow  region  for  the  aspect 
under  consideration. 

(iv)  New  developments  in  the  ray-cptical  method  enable  one  to  treat 
three-dimensional  problems! 20-22 ] ,  and  to  attack  waveguide 
problems  [23-27], 
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To  solve  the  general  scattering  problem  from  a  complex  airborne 
target,  we  follow  a  recent  study  by  Wright  [2].  The  essential  steps  and 
their  appearance  in  this  report  are  described  below: 

(1)  The  target  Is  "resolved"  into  an  ensemble  of  components,  each 
of  which  can  be  geometrically  approximated  by  a  simple  shape. 

This  is  discussed  in  Chapter  2. 

(ii)  The  backscattered  field  from  each  component  is  calculated  by 
the  ray-optics  method.  Details  of  this  step  are  given  in 
Chapter  3. 

(ill)  The  final  step  as  presented  in  Chapter  4  is  to  combine  the 

contributions  from  each  component  in  a  proper  manner,  taking 
into  consideration  the  geometrical  phase  delay,  the  polarization, 
and  the  shadowing  effect.  Ohce  the  total  backscattered  field 
is  ootalned,  the  evaluation  of  various  radar  parameters,  e.g., 
radar  cross  section,  depolarization,  and  glint  errors  [28-30] 
can  be  accomplished  by  straightforward  computations.  Results 
from  those  computations  and  some  comparison  with  experimental 


data  are  given  in  Chapter  5. 


2.  GENERAL  DESCRIPTION 


2.1  The  Targets  and  Coordinate  Systems 

In  this  study  we  will  consider  tvo  specific  targets:  the  BQM-34A 
target  drone  (Figure  2.1)t  and  the  T-33  aircraft  (Figure  2.2).  The 
components  of  the  targets  will  be  modeled  by  perfectly  conducting 
bodies  having  simple  geometries  which  are  described  below: 

(1)  The  ellipsoid  (Figure  2.3),  given  by  the  equation 


2 

x 

“I 

a 


$• 


-  1 


(2.1) 


i 

will  be  used  for  modeling  the  fuselage  and  two  wing  pods. 

(2)  The  elliptical  plate  (Figure  2.4),  given  by  the  equation 


2  2 

■  1  ,  and  s  •  0  (2.2) 

*  b* 


will  be  used  for  modeling  wings,  horizontal  and  vertical  tails. 

(3)  The  open-ended  ,  hollow,  semi-infinite,  elliptical  cylinder 
(Figure  2.6),  given  by  the  equation 

2  2 

~+^«l,  and  z  >  0  (2.3) 

az  o 


4 


3.93m 


Figure  2.1.  BQM-34A  target  drone 


will  be  used  for  modeling  jet  engine  intake  and  exhaust  ducts. 

(4)  The  ogive  (Figure  2.7)  is  generated  by  rotating  an  arc  of  a 
circle  of  radius  R^  about  a  chord  located  a  distance  R^  -  t 
from  the  center  of  the  circle.  The  parameters  (a,  L)  can  be 
calculated  from  the  relations: 


Appropriately,  truncated  ogives  will  be  used  for  simulating 
the  terminations  Inside  the  engine  intake  and  exhaust  ducts. 

To  describe  the  relative  position  of  each  component,  we  will  use 
two  types  of  coordinate  systems:  the  target-centered  coordinate  system, 
with  variables  denoted  by  (x,y,z);  its  origin  is  located  at  the  center 
of  the  ellipsoid  which  models  the  fuselage;  and  its  three  axes  in 
reference  to  the  aircraft  are  oriented  in  the  manner  sketched  in 
Figure  2.1.  In  general,  there  are  as  many  local  coordinate  systems 
as  the  components  that  the  aircraft  are  "resolved"  into.  We  will 
denote  the  local  coordinates  by  (x^n\y^  ,z^) ,  with  n  *  1,  2,  .  .  .  N 
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and  N  being  the  total  nuaber  of  components.  To  describe  a  local  coordinate 


in  reference  to  the  target-centered  coordinate,  the  following  six  parameters 
may  be  used: 

(x<n) ,  ,rt(n)):  Position  of  the  origin  of  the  nth  local  coordi¬ 

nate  (see  left  wing  of  Figure  2.1  as  an  example) 

:  The  angle  represents  the  rotation  of  the 
nth  local  coordinate  about  the  z-axis  of  the  target-centered 
coordinate,  the  angle  6^  about  its  new  y-axls,  and  about 
its  new  x-axls. 

All  rotations  are  defined  in  a  right-handed  sense.  The  explicit 
formulas  for  transforming  variables  from  one  coordinate  system  to 
another  will  be  given  in  Section  4.1. 

In  Tables  2.1  and  2.2,  the  dimensions  and  the  relative  positions 
of  components  that  have  been  used  to  model  the  two  aircraft  under 
consideration  are  presented.  It  should  be  mentioned  that  these 
dimensions  presented  in  Tables  2.1  and  2,2  are  not  exjet.  For  example, 
the  data  for  T-33  aircraft  were  measured  directly  from  a  figure  in 
reference  [33].  Note  that  for  the  BQK-344  target  drone,  there  are 
12  compo'nents  (N  -  12)  in  total.  Obviously,  the  detailed  structure  of 
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OF  THE  BQM-34A 
»  In  meters) 


an  aircraft  Is  far  more  complex  than  that.  The  number  N  ■  12  is  chosen 
merely  as  a  compromise  between  the  accuracy  of  the  modeling  and  the  labor 
of  computations.  Our  computer  program  is  written  in  a  such  manner  that 
if  more  accurate  results  are  desired,  several  additional  components  may 
be  added  conveniently. 

2.2  Scattering  Parameters 

One  of  the  most  important  scattering  parameters  in  a  radar  system 
is  the  backscattered  radar  cross  section  (RCS).  For  an  incident  time- 
harmonic  plane  wave  with  polarization 

t(i)  -  exp  (“ik63  •  r  -  iwt)  (2.6) 

where  (fe^,^,^)  form  a  right-hand  orthonormal  system.  The  backscattered 
field  from  a  complex  target,  in  general,  can  be  written  in  the  following  form 

E(bs)  -  (61E^bs)  +  62E^b8>)  exp  (-ik§3  *  r  -  iwt)  .  (2.7) 

E^bs^  represents  the  amplitude  of  the  reflected  wave  with  same  polariza¬ 
tion,  while  E2ba^  is  the  amplitude  of  the  reflected  wave  with  cross 
polarization.  Finally,  the  RCS  is  defined  as 
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(2.8) 


2  !^b,)l2 

*  |e<»|2 

wher'j  R  Is  the  radar  range.  Another  definition  used  frequently  la  the 
RCS  for  the  cross  polarization, 


o  ■  4irR‘ 

x 


2  iE2 


(bs),2 


I* 


(1)|2 


(2.9) 


The  algebraic  sum  of  a  and  is  sometimes  defined  aa  the  total  RCS  of 


i  the  target,  which  is  a  measurement  of  the  total  reflected  power  with 

reference  to  the  incident  power. 

In  addition  to  RCS,  the  other  parameters  of  interest  in  the  present 
study  are  those  describing  the  target-induced  error  in  the  radar  system. 
While  extensive  discussions  on  the  radar  error  can  be  found  in  the 
literature  (e.g.,  Wright),  we  will  present  here  only  two  definitions  for 
the  orthogonal  errors  caused  by  the  distortion  of  the  reflected  wave  front: 


e  "  Pr  * 


(2.10) 


where  e.  and  e.  are  known  as  the  elevation  and  azimuth  errors,  respectively. 

0  <P 


(?  .P-,P.)  are  the  spherical  components  of  the  time-average  power  flow  p 
r  8  9 
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'>  /•'<"  yt  ~*»*r  »  -r'.T'rw*' 


of  the  backscattered  field 

p  -  Re(j£(bs)  x  S(ba))  (2.11) 

where  Re  is  the  real  operator.  In  Chapter  5  the  parameters  o,  eQ,  and 

e.  will  be  computed  for  each  set  of  input  conditions. 

<? 
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3.  SCATTERED  FIELD  FROM  INDIVIDUAL  COMPONENTS 


There  are  four  basic  geometries,  namely,  ellipsoid,  elliptical  plate, 
open-end  elliptical  cylinder,  and  ogive  that  have  been  used  for  modeling 
the  aircraft  of  Interest.  In  this  section,  the  backscattered.  field 
from  each  of  the  geometries  will  be  determined  Individually  by  the 
ray-optics  method.  In  each  case,  the  problem  will  be  solved  in  its 
respective  local  coordinates  (x^  ,y^  ,z^) . 

Let  us  denote  the  incident  wave  by  (with  time  factor  exp(-iut)  supressed) 

E(i)  -  (0En  +  $E.)  exp(ilc  •  r)  (3.1) 

where  k  ■  (k,0,<fr)  in  the  spherical  coordinates,  and  (0,$)  are  the 
aspect  angles.  The  backscattered  field  can  be  written  in  the  form 

!(ba)  -  (0S0  +  is#>  eXp-^\‘  -k'  -1  ^  .  (3.2) 

Then  the  problem  is  to  determine  (S0,S^)  for  each  geometry  mentioned 
above.  This  will  be  carried  out  below. 

3.1  The  Ellipsoid 

At  high  frequency,  the  main  contribution  from  an  ellipsoid  (Figure  2.3) 
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to  the  backseat ter td  field  comes  from  the  specular  reflection  at  the 


point 


xc  ■  —  sin  0  cos  $ 
a  p 


y„  »  —  sin  6  sin  A 
■'S  p 


z_  ■  — ■  cos  0 
S  p 


(3.3) 


where 


/2 2  2  2  2  2  2  2 

p  “  /a  ain  0  cos  $  +  b  sin  0  sin  $  +  c  cos  0  , 


which  is  the  so-called  "scattering  center"  for  the  ellipsoid.  In  the 


discussion  of  the  shadowing  effect  in  Section  4.2,  the  position  of  the 


scattering  center  relative  to  other  components  of  the  aircraft  determines 


the  posslblity  of  a  contribution  to  the  total  scattered  field  from  the 


particular  ellipsoid  under  consideration.  Thus,  the  Information  in 


Equatiou  (3.3)  Trill  be  used  later.  The  scattering  amplitudes  (S^,S^) 


have  been  worked  out  and  the  results  are 
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(3.4) 


se"£Ke(v?  «P(-i2k^> 

V'Kfifl  $ «P(-12^) 

where 

D  -  [sin2  0  cos2  «  +  (b/a)2  sin2  0  sin2  4  +  (c/a)2  cos2  0]1/2 

(3.5) 

Since  the  ellipsoid  is  a  smooth  body,  there  is  no  cross  polarization 
generated  in  the  scattering  process,  and  furthermore,  the  backscattered 
field  is  Independent  of  the  polarization  of  the  incident  wave. 

3.2  The  Elliptical  Plate 

The  geometry  of  the  elliptical  plate  is  shown  in  Figure  3.1.  If 
the  minimum  radius  of  curvature  of  the  edge  is  large  compared  to  the 
wavelength,  l.e., 

k  b2/a  »  1  ,  (3.6) 

then  the  far- field  coefficients  S.  and  S„  for  vertical  (sf1^  »  0)  and 

$  u  Q 

horizontal  ( E »  0)  polarizations  can  be  asymptotically  expanded  in  a 
9 

-1/2 

power  series  of  k  .  The  domiuant  terms  of  these  series  come  from 
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optical  rays  which  are  backscattered  from  the  scattering  centers 

and  pS2(xs2*?S2'^  on  t^ie  eclRe  of  the  plate.  The  coordinates  of  these  two 

points  are  given  by 


“si  •  -*S2  ■  A  “*  ♦ 


ySl  ’  -yS2  ’  S  sln  * 


ZS1  ■  ZS2  ‘  0 


(3.7) 


where 

A  -  [cos2  *  +  (b/a)2  sin2  *)1/2  .  (3.8) 

According  to  the  geometrical  theory  of  diffraction,  the  first-order 
scattered  field  at  either  of  these  points  may  be  written  in  the  form 


2<bs) 


ei*/4 

fiirk 


TAe1  exp(ikr  +  i¥) 


(3.9) 


where  r  **  the  divergence  factor 

e  ^incident  electrical  field  in  terms  of  the  base  vectors 
T,  B,  and  N  (see  Appendix  A.l) 

¥  **  the  phase  delay. 
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The  divergence  factors  F,  and  r2  associated  with  the  rays  backscattered 
from  Pgl  and  Pg2  can  be  calculated  by  using  Equations  (A, 3.1)  and  (A,  3. 2) 


A  A  A 


in  the  Appendix,  with  p*  ■  — -  ,  0  -  ir/2  ,  d  •*  -k  »  x  sin  0  cos  $ 

aAJ 


+  y  sin  6  sin  <p  +  z  cos  6.  The  results  are 


£  A_3/2(2a  sin  e)”1^2 


(3.10) 


The  column  vectors  for  the  incident  wave  (Ai)pg^>  and  (e*)^^  are 

AAA 

expressed  in  terms  of  the  local  base  vectors.  T^,  at  ,  and 

A  A  A 

T2,  N2,  B2  at  Pg2,  respectively.  The  local  base  vectors  are  illustrated 
in  Figure  3.1,  and  are  defined  by 


AAA 


T^  -  -T2  ■  -x  sin  $  +  y  cos  <► 


■  -N2  -  x  cos  <j>  +  y  sin  $ 


A 


A  A  A 

-  B2  -  -z 


For  the  incident  field  given  in  Equation  (3.2),  we  have 


(®  >PS1  *  E0  «°8  9 

\e0  sin  0, 


(3.11) 
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The  diffraction  matrices  A^  and  A^  for  Pg^  and  Pg2  are  given  by  Equation 
(A. 2. 2) ,  (A.2.5.)»  and  (A. 2./)  in  the  Appendix,  with 

e  -  W/2  , 

a  *  6*  «  8  for  A^  , 

and 

a  *■  0'  «  -0  for  A2  .  (3.13) 

The  diffraction  angles  a,  0'  for  Pgl  and  Pg2  are  illustrated  in  Figure 
3.2.  Thus,  the  diffraction  matrices  are 


(3.14) 
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1 


A2  "  2 


-1  + 


sin  6 


\ 


\ 


1  +  co#2  e 


- 1  + 


sin  8 


sin  8  cos  6 


-  1  + 


sin  6 


sin  6  cos  6  11  f 


dre  8ln2  6  / 


(3.15) 


The  phase  delay  T*  and  T*  for  scattering  fields  at  Pgl  and  Pg2  are 


r..n 


2kA  sin  6 


(3.16) 


Substituting  Equations  (3.10),  (3.12),  (3.14),  (3.15),  and  (3.16)  Into 


A  A 


Equation  (3.9),  we  obtain  in  terms  of  the  base  vectors  T^,  N^, 

Bl*  T2*  N2*  and  B2*  Wrltin*  in  the  £oro  of  Equation  (3.2),  it 


follows  that 


H|  ka____\A-3/2 

ll+  1  ) 

U  ** 8in  el 

J1  +  sinlj 

exp(12kaA  sin  0-1  w/4) 


+  jl  +  alp-g)  exp(-12kaA  sin  0  +  1  ir/4)|  , 


(0  I*  0) 

(3.17) 


The  above  result  does  not  apply  when  0  *  0  and,  therefore,  it  is  not  valid 
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when  the  direction  of  the  incident  wave  is  close  to  the  z-axls.  From 


a  geometrical  optics  viewpoint ,  8-0  defines  a  caustic  where  the 


diffracted  rays  from  all  points  on  the  rim  of  the  plate  converge.  To 


obtain  the  field  in  the  caustic,  we  have  to  use  techniques  other  than 


the  ray-optics  method. 


For  the  exact  case  of  0  -  0,  the  present  problem  can  be  solved 


easily  by  using  the  physical  optics  method  with  the  result 


S<)>  - 

0  1-E0I 


j  ab  k2  ,  (6-0)  . 


(3.18) 


Since  the  first-order  field  obtained  above  is  nonuniform  in  6, 


we  now  introduce  Interpolating  functions  which  yield  results  (3.17)  and 


(3. 18) as  particular  cases.  A  correct  Interpolation  could  be  obtained 


by  asymptotically  expanding  the  exact  solution  of  the  scattering 


problem,  if  such  a  solution  were  available.  A  discussion  jimilar  to 


that  performed  by  Keller  for  the  circular  edge  would  yield  an  axial 


correction  factor  containing  a  combination  of  Mathleu  functions. 


For  numerical  purposes,  it  is  preferable  to  Introduce  matching 


functions  which  are  more  easily  calculated,  such  as  Bessel  functions. 


m 
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We  do  rot  pretend  that  the  formula  thus  produced  is  the  correct  interpo¬ 


lation,  but  onlv  that  it  is  sufficiently  accurate  for  practical  purposes. 

k  comparison  between  Equations  (3.17),  (3.18)  and  the  corresponding 
results  /for  a  circular  disk,  as  well  as  an  inspection  of  the  interpola¬ 
tion  formula  used  for  the  circular  disk,  leads  us  to  suggest  the  following 
expression 


J1(2kaA  sin  8) 


sin  8 


-  IJ^^kaA  sin 


9)]  . 


(3.19) 


The.  above  equation  reduces  to  the  result  for  the  circular  disk  when 

a  ■  b,  to  the  physical  result  Equation  (3.18)  when  8*0,  and  to 

-1/2 

Equation  (3.17)  plus  terms  0[(ka)  ]  when  the  argument  of  the  Bessel 

functions  is  large  with  respect  to  a  unit, 

The  result  in  Equation  (3.17)  is  only  the  first-order  solution  for 
large  ka.  The  higher-order  terras  were  discussed  in  a  separate  report, 
but  they  will  not  be  used  in  the  present  computations. 

3.3  The  Semi-Infinite  Elliptical  Cylinder 

First  let  us  consider  the  case  when  the  cylinder  is  of  circular 


shape,  as  shown  in  Figure  3.3.  For  the  incident  wave  given  in  Equation  (3.1), 


-‘W/vj't'j.uiyi*  r-*+-r*~<~K-~+'~-JKr-->*z^«r-****-pwtrvi.1  «<e3ntWMUnR!CWn 


Z 


Figure  3.3.  The  semi-infinite  circular  cylinder 


the  exact  solution  can  be  obtained  using  the  Wiener-Hopf  technique  [31*32]. 
Omitting  the  details*  the  final  results  are 


E.ka 

6 


l 


eittl *  J2(ka  sin  0) 
n 


(sin  4)  2  n— “  L2(ka  cos  0) 

t  n 


4  ! 
sin  2 


-COS  0 


A2  1 
1  -  A_  J 


ein<^  J  (ka  sin  0)  J'(ka  sin  0)  A 


lE.2ka 

9  _ 

.  -  J*  m  L  (ka  cos  0)  M  (ke  sin  0) 

sin  0  n*-«  n  n 


1  -  A 


(3.20) 


E 


S<J>  “ 


ka  »  ein^  J'(ka  sin  0)  J’(ka  sin  0)  [  4  \  A2  1 

_____  r  n  n  [cos  l  n  I 

.  ®]2  J*  „  M(ka  cos  6)  M  (ka  cos  0)  Icos  6  ~  .  .2  1 

s  ?l  n«-«»  n  n  l  -  A  J 


H !) 


1E.  2ka  »  ein^  J’(ka  sin  0)  J  (ka  sin  0)  A 
0  r  n  n  n 

+  v"  l 


sin  0 


M  (ka  cos  0)  L  (ka  cos  0) 
n  n 


1  -  A_ 


(3.21) 


In  the  above  equations*  the  following  notation  has  been  used 


n  Ln(ka) 

*n  -  iferar 

n 


(3.22) 


The  functions  L  (u)  end  M  (u)  are  analytic  in  the  upper  half  of  the  complex 
n  n 

u-plane  and  are  defined  by 


L  (w)  L  (~u>)  -  itX  J  (1)  H(1)  (X) 
n  n  n  n 

M  (o>)  M  (-w)  -  irX  J’(X)  H^1)’(X)  (3.23) 

n  n  n  n 
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where 


ImX  >  0  for  Imk  >  0  . 


The  above  results  are  too  complicated  to  be  useful.  (Note  that  Ln(u>), 
for  example,  is  explicitly  given  by  an  infinite  integral [ 30] ) .  For  high 
frequencies,  we  may  use  the  approximations 


Ln(w),  Mn(w)  -*■  1  ,  as  ka  ®  .  (3.24) 

Then  the  results  in  Equations  (3.20)  and  (3.21)  can  be  considerably 
simplified 

2  (■§)  J0(2ka  sin  0^  +  sin  0  C0S|2  (7) 

(3.25) 

.g  ka 

%  -  drr  co“2 1!)  Ji<2k* ein  8>  •  <3-26> 


J1(2ka  sin  0) 
2ka  sin  6 


The  above  simple  formulas  are  valid  for  all  6  except  for  6  s  tt/2. 

The  plane  8  -  tt/2  defines  the  shadow  boundary  for  the  specularly 
reflected  wave  from  the  cylinder,  and  the  expressions  for  (SQ ,  S^)  as 
given  in  (3.20)  and  (3.21)  no  longer  hold  since  they  were  derived  from 
simple  saddle-point  Integration.  Fortunately  In  the  present  application, 
we  are  only  interested  in  the  case  with  6  -  0  (head-on  incidence)  and  its 
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immediate  neighborhood  ,i.e.»  8  <  ir/4.  Thus,  this  shortcoming  of  the  formulas 
at  8  -  tt/2  has  no  effect  In  the  later  computations. 

Because  of  Its  application  to  the  radar  return  calculation  from 
an  engine  Intake,  the  problem  of  scattering  by  a  semi-infinite  cylinder 
has  been  studied  extensively  In  the  past  fifteen  years.  For  the  first 
time,  a  simple  analytical  formula  has  been  derived,  l.e.,  Equations  (3.25) 
and  (3.26).  Bowmain  in  1970  gave  a  similar  result  except  that  his 
expression  for  SQ  contains  only  the  first  form  of  (3.25)  and  therefore 
is  not  valid  for  8  -  0  and  its  neighborhood.  For  ka  sin  8  »  1,  the 
results  in  (3.25)  and  (3.26)  are  reduced  to 

1 

s*~  -  iQ  Irsiri  I1 1  sb)  co»<2  -  •'*>  <3-27> 

which  agree  with  the  ray-optics  expression  derived  by  Ross  in  1967. 

As  is  well-known,  the  ray-optics  expression  is  not  valid  in  the  caustic 
region  (8  3  0)  which  is  precisely  the  case  of  interest  in  the  present 
study. 

The  cylinder  considered  above  has  a  circular  shape.  By  proper 


scaling,  an  approximate  formula  for  an  elliptical  cylinder  (with  a  major 


axis  a  and  minor  axis  b  as  shown  in  Figure  2.5),  can  be  derived.  The 
scaling  factor  we  have  selected  is  the  one  between  a  circular  plate 
and  an  elliptical  plate;  it  has  been  rigorously  determined  [3,4].  In 
terms  of.  the  backscattered  field,  the  use  of  the  scaling  factor  leads 
to  this  result: 


.  1/2 


£(bs) (elliptical)  -  (&)  £<ba) (circular  with  radius  aA) 


(3.28) 


where  (aA)  can  be  considered  as  the  "effective  radius"  for  an  elliptical 
cylinder  with 


A  «  >4o82  cf>  +  (b/a)2 sin'*  \  (3.29) 

and  <t>  is  the  azimuthal  angle  of  the  incident  plane  wave. 

3.4  The  Ogive 

The  geometry  of  an  ogive  is  shown  in  Figure  2.6.  The 
scattered  field  from  the  ogive  can  be  obtained  by  applying  the  opt'.cs 
approach  [1].  For  the  present  problem,  we  are  only  interested  in  the 
region  (90°  -  o)  <  6  <  90°,  where  the  scattering  amplitudes  are  given 
by 


S$ 

9 


1  - 


-  t 
sin  9 ) 


1/2 


(3.30) 
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4.  TOTAL  SCATTERED  FIELD 


4.1  Transformations  of  Coordinates 

In  the  previous  chapter,  the  backscattered  fields  from  individual 
components  were  found  and  the  final  results  were  expressed  in  terms 
of  their  respective  local  coordinates.  To  combine  them  properly  for  the 
total  scattered  field,  those  results  have  to  be  transformed  so  that 
they  have  a  common  reference,  namely  the  target-centered  coordinate 
system.  Such  a  transformation  is  straightforward  but  quite  tedious, 
and  a  small  mistake  may  crucially  affect  the  final  result.  Therefore, 
we  will  discuss  various  transformations  in  explicit  detail.  Our  choice 
of  coordinate  systems  is  essentially  the  same  as  that  used  by  Wright 
except  for  some  small  variations. 

For  the  target-centered  coordinate  system  (TCS),  the  x  axis  is 
along  the  longitudinal  direction  of  the  aircraft  with  the  positive 
direction  pointing  from  the  tall  to  the  bead  of  the  aircraft  (see 
Figure  2.1).  The  positive  y  axis  is  directed  to  the  left  side  of  the 
fuselage,  and  the  s  axis  is  along  the  vertical  direction  of  the  aircraft. 

The  other  coordinate  systems  of  interest  are  the  local  coordinate 


* 


w"  »»*"»"'  --roirw*.. 


$ 


M 

I 

r 


// 
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systems  (LCS)  which  are  assigned  to  each  component  of  the  target.  The 
left  wing  is  taken  as  the  example  which  illustrates  the  LCS.  The  left 
wing  is  approximated  by  a  half  of  an  elliptical  plate.  The  origin  of 
the  coordinate  system  is  located  at  the  center  of  symmetry  of  the  ellipse, 
and  the  x  and  y  axes  are  the  axes  of  symmetry  in  which  the  y  axis  is  the 
major  axis  pointing  to  the  left  side  of  the  target  (See  Figure  4.1). 

The  LCS's  for  the  rest  of  the  components  are  assigned  in  a  similar  manner. 

The  LCS  is  obtained  by  making  a  translation  and  three  successive  rotations 
from  TCS.  Let  xt,  yt,  and  zt  be  the  coordinates  of  the  origin  of  the 
local  coordinate  system  referred  to  target-centered  system,  a  be  the  rota¬ 
tion  angle  of  the  LCS  about  its  z  axis,  3  be  the  rotation  angle  about  its 
new  y  axis,  and  y  be  the  rotation  angle  about  its  new  x  axis  (See  Figure  4.1). 
All  three  notations  are  defined  in  a  right-hand  sense.  Let  x,  y,  z  be 
the  coordinates  of  a  point  in  the  TCS,  and  x’,  y',  z*  be  the  coordinates 
of  this  point  in  the  local  coordinate  system.  Then  the  coordinates  of 
the  point  are  related  by  the  following  transformations: 


n 

C12 

ft 

to 

*22 

*32 

33 


(4.1) 


X 


O' ‘origin  of  ioco! 
coordinate  system 

0  ‘origin  of  target 
centered  system 


Figur*  4.1.  Local  coordinate*  for  the  left  wing 
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where 

tn  ■  cos  a  cos  6 

ti2  "  sin  a  cos  3 

t13  -  -  sin  S 

t2i  ■  -  sin  a  cos  7  +  cos  0  sin  3  sin  y 

tj2  ■  cos  a  cos  y  +  sin  a  sin  3  sin  y 

t23  •  cos  3  sin  y 

t^i  ■  sin  a  sin  y  +  cos  a  sin  3  cos  y 

t^2  ■  ~  cos  a  sin  y  +  sin  a  sin  3  cos  y 

t33  »  cos  3  cos  y  .  (4.3) 

In  the  above  transformation,  the  three  successive  rotations  a,  3, 
and  y  have  to  be  made  in  the  order  as  stated;  that  is,  a  first,  3  second, 
and  y  third.  The  corresponding  values  of  xt,  yfc,  zfc,  a,  3,  and  y  for 
components  of  our  targets  are  listed  in  Tables  2.1  and  2.2. 
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Any  point  In  spaca  can  ba  given  in  spherical  coordinates  by  (R,0,$) 


with  respect  to  some  coordinate  system.  The  angles  8  and  4>  which  are 


the  spherical  coordinates  of  the  radar  in  the  TCS  are  usually  referred 


to  as  the  aspect  angles.  The  aspect  angles  are  shown  in  Figure  4.2. 


The  scattered  field  for  the  nth  component  of  the  targets  can  be 
calculated  when  the  aspect  angles  6  ^ for  the  component  in  the 


LCS  are  known.  For  the  given  aspect  angles  8  and  $  in  TCS,  the  aspect 
angles  6^n\  <t>^  in  the  nth  LCS  can  be  calculated  as  follows 


yl 

♦  ^  ■  arctan  ~r 
o 


(4.4) 


and 


,(n) 


arctan 


/*;>2  +  <y;>2 


(4.5) 


where 


s' 


C11  *12  *13 


*21  *22  *23 


*31  *32  fc33| 


cos  t  sin  8 


sin  $  sin  8 


cos  8 


(4.6) 


The  transformation  in  Equation  (4.6)  is  valid  when 


R  »  (*[  +  yj  +  t2)1/2 


(4.7) 
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where  R  le  the  rent*  of  the  reder  referred  to  the  TCS. 

In  calculating  the  beckecetter  A  field  from  eech  component ,  It  Is 
necessary  to  transform  tfc-.,  Incident  field  given  in  the  TCS  to  that  In 
the  LCS,  and  also  to  transform  the  backseat tered  field  in  the  reverse 
direction. 

For  an  Incident  plane  wave  In  Equation  (3.1) ,  the  electric  field 
amplitude  la  expressed  by  a  column  vector 


t1 


(4.8) 


which  Is  written  in  terms  of  base  vectors  In  spherical  coordinates. 
In  the  nth  LCS,  It  becomes 


jl(n)  »  o  T  s 


(4.9) 


where 
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sin  6  cos  $ 
sin  6  sin  t 


cos  8  cos  $ 
cos  ♦  sin  t 
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cos  4 
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with  tjj's  as  in  Equation  (4 

.3)  and 
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*(n) 
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0 
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(4.12) 


Note  that  the  translation  of  the  origins  from  the  target  centered  coordinates 

/ 

Is  neglected  here.  Similarly,  for  a  given  t ^8)  *  in  LCS,  the 
(bs) , (n)  ln  TCS  lfl  glven  by 

|(bs),(n)  .  w  g,  T,  |(bs),(n-i)  (4 


where  U',  T' ,  S'  are  transposes  of  matrices  U,  T,  S  Respectively,  and 


(4.14) 
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with 


■  -xfc  cos  0  cos  ♦  -  yfc  cos  6  sin  ♦  -  *t  >1°  6 
•  -*t  sin  ♦  +  yt  cos  $  .  (4.15) 

The  matrix  W  takes  Into  consideration  the  amplitude  shift  of  the  scattered 

% 

field  of  the  nth  component  caused  by  the  origin  shift  of  its  LCS  from 
the  origin  of  the  TCS.  The  phase  shift  of  the  scattered  field  of  tbe 
component  caused  by  the  origin  shift  is  given  by 

■  -2k(xt  sin  8  cos  ♦  +  yfc  sin  0  sin  ♦  +  *t  cos  0)  •  (4.16) 

In  Equations  (4.15)  and  (4.16),  the  point  (*t»yt»*t)  represents  the 
position  of  the  LCS  origin  referred  to  the  TCS,  and  the  angles  6  and 
4  are  the  aspect  angles  in  the  TCS. 

4.2  Shadowing  Effects 

In  the  present  study,  we  used  the  ray-optics  approach  for  the 
determination  of  the  scattered  field  except  for  the  case  of  caustics; 
the  scattered  fields  are  generated  exclusively  from  the  scattering  at 
the  "scattering  centers."  For  a  given  Incident  angle,  a  scattering 
center  of  a  component  may  or  may  not  be  shadowed  by  other  components  of 
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the  aircraft.  If  it  Is  shadowed,  there  Is  no  contribution  from  that 
scattering  center  to  the  total  scattered  field;  otherwise,  its  entire 
contribution  will  be  included.  To  account  for  the  "all-or-nothing" 
shadowing  effect,  we  will  introduce  a  modulation  function  for 

each  component  of  the  aircraft. 


G(n)(0,4) 


1  if  (e,<(>)  is  in  the  visible  region  of  the  component 


0  if  (6,4)  is  in  the  shadow  region 


(A. 17) 


where  6,  4  are  the  aspect  angles  in  the  TCS.  As  an  example,  the 
shadowed  and  visible  regions  of  the  left  wing  of  the  BQM-34A  target  drone 
are  illustrated  in  Figure  A. 3  when  6  -  90°.  For  the  incident  angle  4 

I 

greater  than  -15°  and  less  than  100°,  the  left  wing  is  in  the  visible 
region  [G^ (0,4)  ■  1];  otherwise,  it  is  in  the  shadow  region 
[G^(6,4)  ■  0].  The  shadow  and  visible  regions  of  the  components  of 
our  targets  are  listed  in  Tables  A.l  and  A. 2. 


A1 


TABLE  4.1  SHADOWING  REGIONS  OF  COMPONENTS  OF  THE  BQM-34A 
(All  angles  are  In  degrees) 


TABLE  4.2  SHADOWING  REGIONS  OF  COMPONENTS  OF  THE  T-33 
(All  angles  are  in  degrees) _ 


1)  Select  aspect  angles  6  and  $  in  the  TCS,  and  an  Incident  field 
In  the  form  of  Equation  (3.1). 

2)  Transform  6  and  $  to  local  aspect  angles  0^  and  using 

Equations  (4.4)  and  (4.5). 

3)  Transform  incident  field  to  local  coordinates  using  Equation  (4.9). 

4)  Compute  backscattered  field  using  one  of  the  equations  in 
Chapter  3. 

5)  Take  into  consideration  the  shadowing  effect,  and  multiply  the 
amplitude  of  the  scattered  field  by  a  proper  modulation  function 
using  Equation  (4.17). 

6)  Transform  the  backscattered  field  to  the  TCS  using  Equation  (4.13). 

7)  Repeat  steps  2)  through  6)  for  all  components. 

8)  Compute  total  scattered  field  in  TCS  by  using 

|0>S>  .  I  g<BS>.<a>  .  „„[„<»>] 

A«1 


where  is  the  phase  shift  in  Equation  (4.16), 

9)  Compute  scattering  parameters  using  Equations  (2.8)  and  (2.9). 
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4.4  Digital  Coaputer  Madding  of  Surfaces 

From  the  above  discussion,  it  is  clear  that  the  two  key  steps  in 
using  the  present  ray-optics  approach  are  (i)  locating  of  scattering 
centers,  and  (li)  determining  which,  if  any,  of  the  scattering  centers 
are  shadowed  by  other  components  of  the  target.  For  simple  geometrical 
bodies,  such  as  an  ellipsoid,  plate,  cylinder,  etc.,  this  task  is  not 

too  difficult.  For  more  complex  geometrical  bodies  that  may  be  con- 

% 

sidered  later,  it  does  not  seem  to  be  analytically  manageable.  In  the 
past  contract  period,  we  have  investigated  a  computer-aided  modeling 
technique,  which  promises  to  generate  arbitrary  surfaces  and  determine 
scattering  centers  and  shadowing  effects  in  an  automatic  fashion.  Some 
of  the  preliminary  results  are  given  below. 

A.  General  Approach  to  Modeling  a  Thrae-Space  Surface 
Two  basic  methodo  exist  for  digital  computer  modeling  of  surfaces. 
The  surface  may  be  represented  by  seme  analytical  function  describing 
the  entire  locus,  or  it  may  be  represented  by  a  data  set  describing 
Individually  sampled  points  on  the  surface.  In  the  course  of  several 
years  of  research  into  digital  display  of  three-dlmensione!  solids,  the 
latter  method  has  been  chosen  as  the  more  flexible  and  efficient. 
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The  general  modeling-  techniques  employed  In  this  work  approximate 
a  three-space  surface  by  a  data  set  consisting  of  an  ensemble  of  planar 
polygons  whose  vertices  lie  on  the  surface  of  interest.  These  polygons 
may  have  any  number  of  vertices  but  must,  however,  be  convex.  For  this 
reason  triangles  or  quadrilaterals  are  most  commonly  used. 

Each  polygon  is  represented  by  a  data  set  which  includes  the 
triplets  of  the  vertices,  a  reflectivity  coefficient,  and  the  three 
coefficients  of  the  general  farm  of  the  plane  in  which  the  polygon 
lies.  These  three  coefficients  are  equivalent  to  the  magnitudes  of 
the  components  of  a  unit  vector  normal  to  the  polygon. 

Once  a  surface  has  been  modeled  by  a  set  of  polygonal  systems, 
programs  are  available  to  display  the  surface  in  a  variety  of  modes. 

The  simplest  form  of  display  consists  of  a  "wire  frame"  or  line  drawing 
of  the  surface.  The  surface  is  represented  by  displaying  lines  along 
all  the  edges  of  the  polygons  in  Figure  (4. A).  Algorithms  are  also 
available  for  producing  a  grey-scale  solid  rendering  of  the  surface  with 
all  hidden  lines  and  planes  removed  [Figure  (4.5)].  These  algorithms 
determine  all  plane  intersections,  all  hidden  or  shadowed  planes,  and 
the  radiant  flux  falling  on  each  plane  because  of  a  movable  radiation  source 


These  algorithm*  essentially  model  the  optical  propertlee  of  the  surface  in 
question  and  then  perform  a  perspective  transformation  for  display  purposes. 

Several  methods  are  available  for  constructing  the  polygonal  data  set. 
Measurements  may  be  taken  directly  from  the  surface  to  he  modeled,  or 
data  points  may  be  generated  from  an  analytical  expression  for  the  surface. 

B.  Polygonal  Modeling  of  a  Surface  Described  by  a  Function  of  Two 
Variables 

» 

A  software  surface  modeling ;  system  has  been  developed  for  operation 
on  the  CDC  1604  computing  system  of  the  Coordinated  Science  Laboratory 
(CSL),  University  of  Illinois.  Parts  of  this  system  will  assemble  a 
polygonal  data  set  from  arrays  of  computed  triplets  lying  on  a  given 
locus.  The  triplets  are  computed  from  the  analytical  expression  of  the 
surface  in  question.  The  values  are  packed  in  three  two-dimensional 
arrays;  for  X,  Y,  and  Z.  The  points  are  related  positionally.  Each 
possible  set  of  four  adjacent  points  is  examined  to  generate  plane  data; 
the  (i,J),  (1  +  1,J),  <l,j  +  1),  «nd  the  (i  +  l,j  +  1)  elements  of  the 
array.  For  every  four-point  set,  two  pairs  of  triangular  planes  are 
possible.  That  pairing  is  chosen  which  produces  the  smaller  dihedral 
angle  between  planes.  This  gives  the  smoothest  approximation  to  the  surface. 
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The  system  then  computes  the  normal  coefficients  and  the  user  may  supply 
a  reflectivity  coefficient.  The  data  set  is  written  on  magnetic  tape  in 
a  format  compatible  with  the  display  algorithms.  This  data  set  may  be 
inputted  to  the  display  systems  or  used  for  further  calculations. 

This  method  gives  a  very  good  approximation  to  the  surface. 

However,  a  few  precautions  are  necessary.  First,  the  function  must  be 
adequately  sampled.  Sampling  frequency  is  determined  by  resolution 
requirements.  The  finer  the  sampling,  the  better  the  fit  of  the  polygon 
function  to  the  surface.  Second,  the  sampling  should  be  finer  at  points 
of  greatest  curvature.  The  polygons  should  not  be  of  equal  area,  rather 
the  difference  in  normal  vector  orientation  between  two  polygons  should 
be  constant,  that  is,  the  angle  between  any  two  polygons  should  be  the 

same. 

For  many  functions  it  is  difficult  if  not  impossible  to  analytically 
derive  expressions  for  curvature.  Even  when  possible,  this  technique  may 
not  be  desirable  because  of  excessive  computation  time.  A  general  procedure 
for  sampling  a  function  so  as  to  produce  a  polygonal  model  reasonably 
sensitive  to  changes  in  curvature  is  to  simply  sample  the  function  in 
a  regular  rectangular  grid  in  the  x,  y  plane.  This  has  the  effect  of 
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producing  larger  planes  as  the  angle  between  the  x,  y  plane  and  the  tangent 
plane  Increases.  This  makes  the  polygon  area  approximately  proportional 
to  the  radius  of  curvature  In  the  region  of  the  polygon. 

A  computer  program  has  been  developed  to  simulate  aircraft  surfaces 
by  an  arbitrary  number  of  ellipsoidal  surfaces.  Each  ellipsoid  is 
specified  by  some  function  z  ■  f (x,y) .  The  coordinate  system  is  centered 
on  the  aircraft  model.  It  is  assumed  that  the  function  contains  terms 
for  rotation  and  translation  which  position  the  ellipsoids  correctly 
with  respect  to  each  other.  The  program  developed  has  been  used  for 
extensive  modeling  of  an  aircraft  composed  of  8  ellipsoids.  The  results 


are  shown  in  Figures  (4.6)  and  (4.7) 


Figure  4.4.  A  "wire  frame"  display  of  an  ellipsoid 


Figure  4.6.  An  aircraft  modeled  by  using  six  ellipsoids  (top  view) 


Figure  4.7.  An  aircraft  modeled  by  using  eight  ellipsoids  (top  view) 


5.  NUMERICAL  RESULTS 

The  numerical  computations  have  been  made  £or  the  two  aircraft: 

T-33,  and  BQM-34A.  In  each  case,  the  RCS,  azimuth  error,  and 
elevation  error  are  calculated  as  a  function  of  the  azimuth  angle  $ 
at  an  Increment  of  1  degree,  while  other  parameters  such  as  elevation 
angle  6,  frequency,  and  polarization  are  fixed. 

5.1  T-33  Aircraft 

In  reference  [33],  the  RCS  of  a  one-half  scale  T-33  model  was 
measured  by  a  nanosecond  short-pulse  radar.  Its  data  provide  a  check 
for  the  present  computer  program.  In  Figures  5.2  and  5.3,  the  experi¬ 
mental  data  and  the  theoretical  prediction  are  presented  when  the 
frequencies  are  8.8  GHz  and  2.9  GHz,  respectively,  with  vertical  polari¬ 
zation  (E0  -  0)  and  0  -  if/ 2. 

It  should  be  noted  that  the  data  given  in  reference  [33]  were 
measured  from  a  half-scaled  model,  while  the  results  presented  in 
Figures  5.2  and  5.3  have  been  increased  6  db  to  compensate  for  the 
scale  factor.  The  agreement  between  the  experimental  and  the  theoretical 
results  is  reasonably  good  except  in  the  neighborhood  of  4  »  0°  (head-on). 
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and  $  •  180°  (tall-on).  In  the  head-on,  or  tail-on  direct Iona,  the 
main  contribution  to  the  scattering  cones  from  the  engine  duct.  The 
fact  that  good  agreement  is  not  obtained  in  these  two  directions  indi¬ 
cates  the  need  of  an  improved  model  for  the  engine  duct.  This  task  will 
be  carried  out  in  the  next  contract  period. 

In  addition  to  the  case  6  ■  90°  as  discussed  above,  the  experi¬ 
mental  data  for  the  T-33  are  also  available  when  the  aircraft  is  tilted 
at  an  angle  Y  from  the  horizontal  plane  (Figure  5.1).  The  corresponding 
aspect  angles  (8 ,  <p)  due  to  such  a  tilt  can  be  found  from  the  transfor¬ 
mation  in  Equation  (4.1)  with  0  ■  -V,  a  «  0,  and  y  *  0.  For  the  frequency  8.8  GHz 

o  o  • 

with  Y  ■  5  ,  and  for  the  frequency  2.9  GHz  with  Y  »  5  and  Y  *  10  ,  the 

experimental  and  the  theoretical  calculations  are  presented  In 

Figures  5.4,  5.6,  and  5.7  respectively.  The  agreement  is  reasonably  good, 

but  not  as  good  as  that  presented  in  Figures  5. 2  and  5.3.  This  is  due 

to  the  facts  that  the  experimental  model  of  the  T-33  is  more  complex 

in  structure  than  the  model  described  in  Table  2.2,  and  that  the  difference 

In  structures  becomes  more  predominant  when  the  aspect  angle  9  is 

different  from  n/2.  The  accuracy  can  be  improved  if  more  components  are 
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used  in  Table  2.2.  However,  this  can  be  achieved  only  at  the  expense 


of  computational  labor. 


S' 

& 


5.2  BQM-34A  Target  Drone 

This  aircraft  was  also  studied  by  Wright  who  used  ellipsoids  to 

i 

approximate  the  various  components  of  the  aircraft.  In  Figures  5.7,  5.8, 
and  5.9,  we  present  a  comparison  between  his  results  and  ours  for  the 


RCS,  eQ,  and  e.  as  functions  of  the  azimuthal  angle  $  at  frequency  1  GHz, 
with  vertical  polarization  Incidence  (E.  *  0),  and  6  *  ir/2.  The 
agreement  is  good  except  in  the  neighborhood  of  head-on  and  tail-on 


directions,  because  the  scattering  from  the  duct  is  taken  into  consideration 


in  our  computation  (See  Sections  3.3  and  3.4),  while  it  is  not  included 


in  Wright *8  computation. 


To  study  the  frequency  dependence  of  the  scattering  parameters, 


more  computations  were  made  by  changing  the  frequency  from  1  GHz 


(Figures  5.10,  5.11,  and  5.12)  to  5  GHz  (Figures  5.13,  5.14,  and  5.15), 


and  to  8  GHz  (Figures  5.16,  5.17,  and  5.18),  while  the  other  parameters 


are  kept  the  same  as  the  ones  used  in  Figures  5.7,  5.8,  and  5.9.  As 


may  be  expected,  the  RCS  is  not  frequency-sensitive  in  the  broadside 
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direction.  In  that  direction,  the  a* In  contribution  to  the  scattering 
cones  from  the  fuselage,  which  is  modeled  by  an  ellipsoid  in  the  present 
study.  As  given  in  Equation  (3.4),  the  backscattering  from  an  ellipsoid 
is  independent  of  frequency. 
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AZIMUTH  ANGLE  IN  DEGREE 


.2,  RCS  in  db  per  square  meter  of  a  full-scaled  T-33  aircraft 
as  a  function  of  azimuth  angle  $  for  a  fixed  0  ■  90°,  for 
a  vertically  polarized  incident  plane  wave  (E  «  0) ,  and 
at  frequency  8.8  GHz.  The  solid  line  represents  the  5°~ 
average  of  the  theoretical  computation.  The  dotted  line 
represents  the  5°~medlan  of  the  experiment  reported  in 
reference  [33]. 
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Figure  5.3.  RCS  in  db  per  square  meter  of  a  full-scaled  T-33  aircraft 
as  a  function  of  azimuth  angle  $  for  a  fixed  8  ■  90a.,  for 
a  vertical ly  polarized  incident  plane  wave  (E.  *  0) ,  and 
at  frequency  2.9  GHz.  The  solid  line  represents  the  5°- 
average  of  the  theoretical  computation.  The  dotted  line 
represents  the  5°-median  of  the  experiment  reported  in 
reference  [33] . 
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Figure  5. A.  KCS  in  db  par  square  metes  of  a  full-scaled  T-33  aircraft 
as  a  function  of  azimuth  angle  $  for  a  tilt  angle  ¥  *  5°, 
for  a  vertically  polarized  incident  plane  wave  (K0  *  0) , 
and  at  frequency  8.8  GHz.  The  solid  line  represents  the 
5°-average  of  the  theoretical  computation.  The  dotted 
line  represents  the  5°-aedian  of  the  experiment  reported 
in  reference  [33]. 


RZIMUTH  ANGLE  IN  DEGREE 


Figure  5.5.  RCS  in  db  per  square  meter  of  a  full-scaled  T-33  aircraft 
as  a  function  of  azimuth  angle  for  a  tilt  angle  ¥  «  5°, 
for  a  vertically  polarized  incident  plane  wave  »  0) , 
and  at  frequency  2.9  GHz.  The  solid  line  represents  the 
5°-average  of  the  theoretical  computation.  The  dotted 
line  represents  the  5°-medlan  of  the  experiment  reported 
in  reference  [331. 
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Figure  5.6, 


AZIMUTH  ANGLE  IN  DEGREE 


RCS  in  db  per  square  meter  of  a  full-scaled  T-33  aircraft 
as  a  function  of  azimuth  angle  $  for  a  tilt  angle  ¥  -  10° 
for  a  vertically  polarized  incident  plane  wave  (E  »  0), 
and  at  frequency  2.9  GHz.  The  solid  line  represents  the 
5°-average  of  the  theoretical  computation.  The  dotted 
line  represents  the  5°-median  of  the  experiment  reported 
in  reference  [33]. 
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Figure  5.7.  RCS  In  db  per  square  meter  of  the  BQM-34A  drone  as  a 
function  of  azimuth  angle  $  for  a  fixed  9  ■  90°,  for 
a  vertically  polarized  Incident  plane  wave  (E.  »  0) , 
and  at  frequency  1  GHz.  The  solid  line  represents  the 
5°-medlan  computation  of  our  model.  The  dotted  line 
represents  the  5°-medlan  computation  of  Wright's  model 
using  ellipsoids. 
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Figure  5.8.  Azimuth  error  (e.)  in  meters  of  the  BQM-34A  drone  as  a 
function  of  azimuth  angle  $  for  a  fixed  8  -  90°,  for  a 
vertically  polarized  Incident  plane  wave  (Efl  •  0) ,  and 
at  frequency  1  GHz.  The  solid  line  represents  the  5°- 
median  computation  of  our  model.  The  dotted  line  repre¬ 
sents  the  5°-median  computation  of  Wright's  model  using 
ellipsoids. 
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Figure  5.9.  Elevation  error  (e.)  in  meters  of  the  BQM-34A  drone  as  a 
function  of  azimutn  angle  $  for  a  fixed  6  •  90°,  for  a 
vertically  polarized  incident  plane  wave  (Eg  *  0),  and 
at  frequency  1  GHz.  The  solid  line  represents  the  5°- 
median  computation  of  our  model.  The  dotted  line 
represents  the  5°-median  computation  of  Wright's  model 
using  ellipsoids. 
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Figure  5.11.  Azimuth  error  (e, )  in  meters  of  the  BQM-34A  drone  as 
a  function  of  azimuth  angle  $  for  a  fixed  9  -  90°, 
for  a  vertically  polarized  incident  plane  wave  (E@  *  0) , 
and  at  frequency  1  GHz  0 
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.12.  Elevation  error  (e.)  In  meters  of  the  BQM-34A  drone  as  a 
function  of  azimuth  angle  $  for  a  fixed  6  -  90°,  for  a 
vertically  polarized  Incident  plane  wave  (E@  *  0) ,  and 
at  frequency  1  GHz 
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Figure  5.14.  Azimuth  error  (e.)  in  meters  o£  the  BQM-34A  drone  as  a 
function  of  azimuth  angle  $  for  a  fixed  6  *  90°,  for  a 
vertically  polarized  Incident  plane  wave  (E0  ■  0),  and 
at  frequency  5  GHz 
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Figure  5.15.  Elevation  error  (e„)  in  meters  of  the  BQM-34A  drone  as  a 
function  of  azimuth  angle  $  and  a  fixed  6  ■  90°,  for  a 
vertically  polarized  incident  plane  wave  (Eg  -  0) ,  and 
at  frequency  5  GHz 
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Figure  5.17.  Azimuth  error  'e  )  in  meters  of  the  BQM-34A  drone  as  a 
function  of  aziir.th  aagle  $  for  a  fixed  0  ■  90°,  for  a 
vertically  polarized  Incident  plane  wave  (EQ  *  0) ,  and 
at  frequency  8  GHz 
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Figure  5.18.  Elevation  error  (e.)  in  meters  of  the  BQM~34A  drone  as  a 
function  of  azimuth  angle  4  for  a  fixed  8  *  90°,  for  a 
vertically  polarized  incident  plane  wave  (Eg  »  0) „  and 
at  frequency  8  GHz 
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6.  CONCLUSION 


\ 


V 

> 


i. 


In  thia  work,  we  have  derived  a  relatively  simple  and  analytical 
model  for  a  complex  airborne  target.  The  results  for  the  radar  scattering 
parameters,  l.e.,  RCS  and  angular  scintillation,  are  given  in  terms  of 
explicit  simple  functions  and  therefore  can  be  calculated  in  an  efficient 
manner  with  the  aid  of  a  computer.  To  test  the  accuracy  of  our  model, 
the  theoretical  results  are  compared  with  the  experimental  data  for  T-33 
aircraft.  The  agreement  is  generally  good  (within  several  db)  except  in 
the  head-on  and  tail-on  directions.  In  these  two  direction,  our  simple 
model  using  a  semi-infinite  cylinder  with  an  ogive  as  the  engine 
termination  to  approximate  the  engine  duct  is  not  adequate;  and  future 
work  to  develop  a  more  realistic  model  is  needed. 
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APPENDIX 


A.l  Diffraction  Matrix  for  a  Wedge 

Consider  a  perfectly  conducting  wedge  of  half-angle  ft  whose  edge 
coincides  with  the  z  axis  of  a  Cartesian  coordinate  system  (x,  y,  z)  and 
whose  surfaces  are  defined  by  the  equations  y  -  ±  x  tan  ft,  x  >  0, 
and  -<*>  <  Z  <  ®  .  The  geometry  of  the  wedge  is  shown  in  Figure  A.  1.1. 

The  problem  of  a  plane  wave  incident  or  on  a  plane  perpendicular  to  the 
Z  axis  with  either  its  electric  or  magnetic  vector  parallel  to  the  edge 
is  a  classic  one;  the  exact  solution  is  known,  and  from  the  basic  results 
it  Is  a  trivial  matter  to  deduce  the  corresponding  solutions  for  plane 
waves  incident  at  an  angle  it-6  to  the  Z  axis. 

According  to  Keller's  geometrical  theory  of  diffraction,  we  define 

A  A  A  A 

a  set  of  edge  fixed  base  vectors  T,  N,  and  B,  where  T  is  a  unit  vector 
parallel  to  the  edge,  N  is  a  unit  vector  normal  to  the  edge  and  pointing 

A 

away  from  the  wedge,  and  B  is  the  unit  vector  binormal  to  the  edge  and 
pointing  into  the  shadowed  half  space.  The  direction  of  T  is  chosen 
such  that  T  *  N  x  B  in  a  right-handed  system.  In  general,  we  choose 
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Figure  A. 1.1.  The  geometry  of  Che  wedge 
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A  J 

The  vectors  e  and  a  ere  related  by 

ed  ■  Ae*  (A. 1.5) 

where 

/-(x  -  y) 

A  ■  -j  (x  -  y)  cot  8  sin  0*  (x  +  y)  cos  0*  cos  a  (x 

\-(x  -  y)  cot  8  cos  0*  (x  +  y)  sin  0'  cos  a  (x 

* 

x  -  —  sin  —  [cos  —  -  cos  —  (a  +  O')!*1 
n  n  n  n 

y  *  ~  sin  ^  [cos  7  4-  cos  —  (tt  +  25  -  a  - 

n  n  n  n 

with 

n  -  Z(1  -  7)  .  (A. 1.8) 

If 

*1 

Both  e  and  e  are  coluan  vectors  with  coaponents  along  the  base  vectors 

A  A  A 

T,  N,  and  B. 

A.  2  Diffraction  Matrix  for  a  Thin  Edge 

The  perfectly  conducting  half  plane  shown  In  Figure  A. 2.1  Is  a 
special  case  of  the  wedge  in  Figure  A. 1.1,  with  half-anglo  ft  equal  to 
zero.  Us  define  the  base  vectors  T,  K,  and  B  as  shown  in  Figure  A. 2.1. 


+  y)  cos  0'  sin  a 

+  y)  sin  G '  sin  a 
(A. 1.6) 

e*)f 1  (a. 1.7) 
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The  Incident  and  diffracted  field  can  be  written  in  the  form  of 


''i  *d 

Equation  (A. 1.2)  and  Equation  (A. 1.3).  The  vectors  e  and  e  are  related 


by 


Ad  aa1 

e  ■  Ae 


(A.  2.1) 


where 


-(-x  +  y) 


A  ■  I  (-a  4-  y)  cot  0  sin  0  ’  (-x-  -  y)  cos  0 '  cos  a  ( -x  -  y)  cos  0  ’  sin  a 


(-x  +  y)  cot  0  cos  0*  (-x  -  y)  sin  0’  cos  a  (-x  -  y)  cos  0'  sin  a 


(A.  2. 2) 


with 


v  „  If 

—  <  a  <  — 

2  -  -  2 


f ;  •’  ir 


(A.  2. 3) 


If  the  observation  point  p  Is  not  near  the  boundary  of  the  shadow  region , 


where  0*  •  a  +  *  ,  l.e. ,  if 


cos 


0'  -  a 


t'2itr  si  a  0  »  1  , 


(A. 2. 4) 
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2  cos 


1 _ 

0*  -  a 


(A.  2. 5) 


If  p  Is  not  near  the  boundary  of  the  reflected  field,  where  0*  **  -a, 


i.e.,  If 


0’  +  a 


/2kr  sin  6  »  1  , 


(A.  2. 6) 


2  sin 


FTa 


(A.2.7) 


A. 3  The  Divergence  Factor  for  the  Diffracted  Ray  from  a  Curved  Edge 


By  applying  the  geometry  theory  of  diffraction,  the  divergence 


factor  for  the  diffracted  field  from  a  curved  edge  can  be  written  as 


L(p  +  s)sj 


(A. 3.1) 


where  p  is  the  principal  radius  of  curvature  of  the  wave  fronts 


diffracted  by  the  edge  and  S  is  the  distance  from  the  edge  to  the 


observation  point.  The  radius  of  curvature  of  the  wavefront  is 
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p 


(A.3.2) 


.zeliisLg. 

A  A  i 

(d  -  i)  •  n 


!  f 


where  p*  is  the  radius  of  curvature  of  the  edge,  8  is  the  angle  between 

A 

the  incident  ray  and  the  tangent  to  the  edge,  d  is  the  unit  vector  in 

A 

the  direction  of  the  diffracted  ray,  and  n  is  the  unit  vector  in  the 
direction  of  the  principal  nonaal  to  the  edge. 

The  principal  norsal  to  an  edge  is  the  normal  lying  in  a  plane 
containing  the  edge  at  the  location  of  Interest. 
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